Diffusion on asymmetric fractal networks 
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We derive a renormalization method to calculate the spectral dimension d of deterministic self- 
similar networks with arbitrary base units and branching constants. The generality of the method 
allows the affect of a multitude of microstructural details to be quantitatively investigated. In 
addition to providing new models for physical networks, the results allow precise tests of theories 
of diffusive transport. For example, the properties of a class of non- recurrent trees {d > 2) with 
asymmetric elements and branching violate the Alexander Orbach scaling law. 



The behavior of random walks on networks has been 
linked to a wide range of interesting phenomena in 
many disciplines [H. In recent decades the study of the 
properties of fractal and, more recently, complex net- 
works 0, 0] has revealed a range of "anomalous " behav- 
ior, characterized by power law scaling with non-integer 
exponents or dimensions. In general, interactivity on net- 
works can be codified in a matrix, which in many appli- 
cations (e.g spring systems, random walks and conduc- 
tivity) takes a Laplacian form [a . Macroscopic proper- 
ties, such as phase-transitions [Jill, and mass and elec- 
tronic transport [1,0] are connected to the eigenspectrum 
of the matrix which is characterized by the spectral di- 
mension d. The dimension d can be found using the 
site-independent [1| asymptotic form of the probability 
Po{t) '-^ i"''/^ that a walker returns to its origin at time 
t. 





Although many networks are random, simple models 
such as deterministic fractals, combs and trees have long 
been of theoretical interest because their properties can 
be analytically studied; the results used to test exist- 
ing results, motivate new theories, mimic physical net- 
works, or simply to provide insights into the nature of 
processes on networks. In this paper we use renormaliza- 
tion ideas to study d for classes of deterministic 
networks that include asymmetric elements and branch- 
ing constants (See Figs. [T] & [2]). The models may there- 
fore be able to mimic asymmetric branching and chirality 
of polymers, as well as aspects of complex networks which 
exhibit high variability in their degree distribution 

In a seminal paper, Alexander and Orbach (AO) 
argued that d was related to the fractal dimension df 
and random walk dimension du, by the simple formula 
d = 2df/dw This result is one of the most important 
relations in the theory of fractals. All analytic results 
to date confirm the law, and Teles has proved that it 
holds for recurrent {d < 2) loopless networks which sat- 
isfy a technical smoothness condition on the electrostatic 
potential which arises in the definition of the network 
resistance. No proof has been given purely in terms of 
geometric and topological characteristics of networks. To 
demonstrate the renormalization method we calculate d, 
df and d^, for a class of networks which violates the AO 
law. In a related study [l3| we explain the discrepancy 
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FIG. 1: Examples of networks. In (a) the base unit is a 
periodic comb with seven teeth and one diagonal strut. In 
(b) the network has branching constant u = 2 and the base 
unit has the shape of an asymmetric staple. 

in terms of anisotropic diffusion on the network. 

The method of iteratively constructing the infinite net- 
work is shown in Fig.[2l The tree is made by taking a base 
unit, doubling its size, and attaching Ui copies of the re- 
scaled unit to each of the two end points of the base and 
so forth. There are a number of generalizations which 
can be readily incorporated, allowing the affect of nu- 
merous topological network characteristics to be probed. 
First, the base unit can be an arbitrary sub- network with 
L left, and R right nodes (The tree example in Fig.[2]has 
L = 1 and R = 2). li L > R, then L — Rot left nodes can 
be fused to form a network with L ^ R. U L < R there 
are 'R choose L' distinct ways of connecting the first it- 
eration to the base; a different branching constant Ui can 
be associated with each possibility. Second, each of the 
elements in the nth iteration of the tree can be replaced 
by the nth generation of a fractal with two specified end 
nodes (e.g. two corners of a Sierpinski triangle). Dif- 
ferent sub-classes of these networks have been studied 
previously [13, El, H [13, [3 ■ 

In standard diffusion the probability that a random 
walker is at a point x at time i on a pipe is governed 
by the equation pt — Dp^x- To formulate the equations 
governing the probability on a network, we first solve the 
diffusion equation on a single bar for arbitrary Dirichlet 
conditions and a homogenous initial condition; p (0, t) = 
Poit), p{b,t) — pi{t) and p(a;,0) = 0. The Laplace 
transform (LT) of the solution is [l9| 
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The flux entering the bar at the node a; = and x = b 
is qo{t) = —Dpx (0, t) and qi{t) — Dpx {b, t) respectively. 



2 




(c) 



(d) 



FIG. 2; (a) The base unit, (b) the second iteration of the net- 
work is obtained by attaching multiple copies of the rescaled 
base unit to the end points of the base unit, (c) a regular pipe 
and (d) an iterated fractal. 

For networks it is useful to express these fluxes in terms 
of the concentrations Ci at either end. This gives rise to 
the "flux-concentration" equations 



tanh(6^s/L») 


jo 




" 1 -a' 




Co 








-a 1 




Cl 



jo 




" 1 -/■ 




Co 


.ji . 




1 . 




. '^1 . 



where a = sec\l{h^J s/ D). For simplicity we set D to 
unity. In general, the bar can be replaced by the genera- 
tor of a fractal (see Fig.[2lld)). The above matrix is then 
a special case of the more general equation 
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Here / is the LT of the first passage time density between 
the ends of the generator and g is the LT of the probabil- 
ity that a walker released at the left end, is at its origin at 
time t. These functions are defined for a generator with 
reflective boundary conditions at either end (l^L Appx. 
A]. For a pipe, f = a and g — {ta,nh.{by/s)y/s)~^ while 

~ -3a2)-i 
I Appx. 



for the generator of a determinstic tree / = ct (4- 
and g = {4- 3a^){A{l - a'^))-Hanh{b^)/y^ ^ 
B]. ^ 

Using these relations, diffusion on an arbitrary network 
can be formulated as a system of algebraic equations [l^ . 
To link two or more pipes, the concentrations at the com- 
mon node are taken to be equal. Mass conservation at a 
node is enforced by requiring that the total flux entering 
all the pipes at that node sum to zero. Using the node 
labels in Fig. E^a), the equations for the example are 



(l-,/').9j 



where c = [co, Ci, C2, C3, 04]^ etc. The number of equa- 
tions for any sub-network with L left end nodes, R right 
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end nodes and K interior nodes can be reduced by set- 
ting jk = (fc = I..K) and solving the under determined 
system to give 
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Here E is an {L + R)x {L + R) matrix and the subscripts 
in the above matrices denote their size (e.g. Jli is an L 
by 1 matrix). 

It is useful to rewrite equation ([1]) as Fj = c, where 
F = (1 — f^)gE^^ is a matrix of Green's functions. The 
equation can be written in terms of sub-matrices as 



(2) 



Tree networks of the type shown in Fig. [T] and [2] can 
be iteratively constructed from these reduced networks. 
Carrying out the reduction process for a network of arbi- 
trary iteration gives a system of L + R^^ equations, where 
i?e is the number of end nodes. If no- flux conditions are 
applied at these nodes (ji = 0, i = 1, . . . , i?e), the equa- 
tions reduce to an L x L system. For the infinite network 
these equations are written as 
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where the hat denotes a quantity associated with the 
infinite network. To completely specify the problem, L 
boundary conditions have to be applied at the left end 
nodes. At the origin, which is defined to be the first of 
the L nodes, an instantaneous unit source qoit) = 5(t) 
is applied (so jo{s) = ^((7o(i)) = !)■ No-flux conditions 
(i.e. ji — 0) are applied to the remaining L — 1 nodes. 
Po(i) is then given by the inverse LT of the first entry of 
the matrix F. If the origin is shifted, its concentration is 
given by the corresponding diagonal entry of F. 

Because of the self similar nature of the tree, the matrix 
F can be studied through renormalization. The first step 
is to consider the equations for the infinite network which 
has been doubled in size. This is achieved by replacing h 
in dSl) by 26 to get 

~ ^g(sech(26v/i)) ~ 

tanh(26Vi) " ^ ^^^^ 

The tilde is used to denote quantities associated with the 
rescaled network. Comparison between the two Green's 
functions in Eqns. (O and (g]) reveals that F(s) = 2r(4s). 

The original network is reconstructed by patching the 
rescaled network to the initial finite network with L left 
nodes and R right nodes. If L = i? this is simply a 
matter of taking Jjn +uJli = 0, where u is the branching 
constant (e.g., u = 2 in Fig. [H (b)). For the example in 
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Fig.[2l(-R = 2, -L = 1), the mass conservation condition is 
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In general, the conservation of mass condition can be 
written as 



jRi = -\FRR{n,f-\^s))Cm, 



(6) 



where u = (ui, U2, ■ ■ ■) are the branching constants. The 
form of Frr is dependent on the structural details of the 
connections. In general Ffm depends on the elements of 
individually. 

Now Eq. dS]) is used in the second equation of to 
express Cjn in terms of J^i, which can then be used in 
the first equation of ^ to get F'^ Jli = Cli, where is 
calculated by algebraic methods from the terms above. 
The renormalization is completed by noting that must, 
by construction, be identical to F given in Eq. Q. This 
leads to the closed equation 



The result can be extended to the case where each pipe is 
replaced by a two-ended fractal. Subsequent generations 
of the network are defined by increasing the iteration of 
the fractal rather than doubling the length of the pipes. 
For an infinite network, Eq. ([3]) becomes 
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Increasing the generation of each fractal within the in- 
finite network gives a structure whose equations (analo- 
gous to Eq. (HI)) are 

Jli = (5*(1 - r^))-'EinCLi = f-'CLi- (9) 

Here /* and g* are the first passage time and Green's 
function associated with the two ends of the first itera- 
tion of the generator with reflective boundary conditions 
applied at both exterior nodes. 

The functions /(s) and /*(s) are related by f{p{s)) — 
/*(s) where p(s) is the renormalization of first passage 
time function [TO, [lOl- For example, the deterministic 
tree has p{s) = arcsech(cr^/(4 — 3ct^))^ Applying 
this relation to the two infinite Green's function forms 
given in © and ^ gives F(s) = g*{s)f{p{s))/g{pis)). 
Proceeding as above gives 



g*Fj,],{u,f{p{s))-^)) 



Fll — Fifl I '■ — , 1- Tiijij Trl- 

(10) 

Regarding the pipe as a fractal with p{s) = 4s recovers 
the prior case [Eq. ([7])] as g*{s)/g{p{s)) = 2 and g*{s) = 
{ta,nh{2by/s)^/s)~^ . Note that the renormalization of the 



tree with fractal elements is undertaken with respect to 
time, rather than space. 

The spectral dimension d is found by considering the 
asymptotic behavior of the first element of the matrix F, 
denoted by 7, which has the form 



7(s) = { kos^^^-Hnis) + 
ko + kis^/"^-^ + 



d<2 
d = 2,4,6,... 

d > 2,d ^ 4,6, , 



(11) 



The Green's function of the zeroth and first iteration have 
the asymptotic forms 

g^l/{Ms) + 0{s^) g* ^l/{M*s) + Ois°), (12) 

where M and M* are the mass of the zeroth and first 
iterations of the generator. The mass rescaling factor 
of the fractal between successive iterations is given by 
A = M*/M. For s ^ 0, p{s) = as + O(s^), where a 
is the time rescaling factor. The ratio of the mean first 
passage times between the left and right end nodes of the 
n — 1th and nth iterations of the generator is equal to a. 

As an example, consider the infinite structure in 
Fig. Wijo) ■ The matrix Frr required for the evaluation of 
Eq. (fTO|l is determined (from Eq. ([5])) to be 



Frr — 



1 



7(p(s)) 



ui 

U2 



(13) 



The spectral dimension can be calculated by substituting 
(fT5)) into (|10p and then using the asymptotic forms (fTT|) 
and (fT2l) . There are several distinct cases. liA{ui+U2) — 
a < then d < 2 with d = 21og^(A(ui + U2)), while if 
A{ui + U2) — a > then d> 2 with 



2 log. 



A ((ui + U2) a; H- 6 U2U1Y 



A:Ui^U2 + (16U2^ -\- -\- I2XU2) Ui + X-^U2 I 

where x > is found by solving the quadratic equation 

{{ui + U2) A — a) 
-\-2 X {2iAuiU2 — a(3u2 -I- 2 ui)) — 20 U2Uia. 

Interestingly, the latter formula for rf > 2 reduces to d = 
21ogQ(^(ui + U2)) when ui = 2u2- 

The AO law can be readily checked for a pipe {A — 2 
and a = 4). At large t and r the random walker's behav- 
ior will be dominated by diffusion on long pipes which im- 
plies du, — 2. A simple renormalization argument shows 
that the mass M{r) of the network at a large distance r 
from the origin is given by M{r) — 4: + 2uiM {{r — 3) /2) + 
2u2M{{r — 2)/2). Assuming M{r) ^ r'^f and taking ex- 
pansions gives df = ln(2(ui +U2))/ln(2). For d <2, and 
for the exceptional case of d > 2 (ui = 2^2), the pipe 
networks therefore obey the AO law. Exluding the case 
ui = 2m2 the results show that d ^ 2df/dw if d > 2. 
Computations confirming the analytic values of dm and 
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FIG. 3: Simulation of d and dw The left hand axis gives 
po{t); we calculate d = 4.37 (Best Fit: 4.36). On the right 
hand a:xis we plot \/r^t) against t^^^ confirming du, = 2. 

d for ui = 1 and U2 — 20 are shown in Fig. [31 The 
properties of the tree network with a general fractal gen- 
erator are qualitatively similar; it is apparent that the 
form d — 2log^{A{ui + U2)) is just the AO law. 

Our main example demonstrates that networks with 
asymmetric characteristics can have properties which vi- 
olate the AO law. Interestingly these unusual networks 
provide an example of a system in which local micro- 
scopic details affect macroscopic properties. This is not 
expected in the scaling theory of dynamic phenomena; 
in general the network properties should be determined 
by power law exponents while microscopic details are 
relegated to pre-factors. This idea is demonstrated by 
considering the addition of a bar between nodes 1 and 
3 of the tree shown in Fig. [21 As expected, this does 
not change df or (A and a remain unchanged). In 
contrast, J is changed significantly and AO now holds if 
ui = 4u2 when d > 2. 

We have argued in Ref. [l^l that the AO law holds if 
each site a distance r from the origin is explored approx- 
imately uniformly, and numerically demonstrated that 
this is not the case for d > 2. The fact that AO holds 
for the tree if d < 2 demonstrates that asymmetry alone 
is not a sufficient condition for anisotropic transport for 
recurrent walks {d < 2). The generality of the renor- 
malization method proposed should assist in the identi- 
fication of geometric and topological characteristics as- 
sociated with anisotropic (non AO-law) transport. In 
particular, the algebraic formulation allows a number of 
results to be stated for arbitray sub- networks. For exam- 
ple, it is possible to prove [2l[ that the AO law is obeyed 
for the class of networks considered here if d < 2 or for 
d>2 when L = R = 2. 

In Ref. [13] it was shown that the resistivity exponent 
C of a network satisfies the law C = %-(2 — d), which 
differs from the conventional relationship [3| C, = dw — df 
if d 7^ 2df /dw A non-trivial test of the general formula 



is provided by replacing the pipe of the base unit by a 
deterministic tree generator [A = A, a = 8). This gives 
d = 3.58, d^ = 3 and hence Q — —2.37, which matches 
the exact value 

We have derived a method for finding d for self-similar 
trees based on an arbitrary sub-network with fractal el- 
ements. The ability to incorporate loops, branching and 
asymmetry allows the affects of these characteristics to 
be studied, and increases the classes of physical networks 
that can be quantitatively modeled. An advantage of 
using the continuum diffusion equation (compared with 
discrete-space walkers) is that non-integer bar lengths 
can be studied by taking p{s) = 4^3, with v real and pos- 
itive. This allows d to be continuously varied and careful 
study of the transition between recurrent and transient 
random walks (at d = 2), as well as the transition from 
fractal to inhomogeneous trees {v 0) (e.g. the Bethe 
lattice). 
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